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^ I Abstract 

f^ ' The efficiency of a simple model of crossflow fan is maximized when the geom- 

• , etry depends on a design parameter. The flow field is numerically computed 

^ ' using a Galerkin method for solving a Poisson partial differential equation. 
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C^ ; 1 Crossflow fans 

'nT ' A crossflow, or tangential, fan (fP) is a ventilating structure where air is moved 

■ , from an inlet zone to an outlet one using an impeller that forces the air to follow 

C^^ ' trajectories tangentially to the ideal circumferences concentrical to impeller it- 



self, so that fluid goes across fan on planes that are normal to impeller rotation 
axes. The inlet and the outlet zone are physically separated by a fixed object, 
called vortex wall (see FigH]). 

Crossflow fans have the important property to get relatively high values of 
^^ , flow rate with small geometrical dimensions. Two typical engineering questions 

$H ' are how to design their housing and where to place the vortex wall between the 

inlet and the outlet zone ([5]). In this work, we try to get an answer to the 
second question using a simple two-dimensional geometry in the Mathematica 
numerical environment. 

2 The geometry 

We consider a squared housing D with edge of length L. The impeller is ideally 
represented by a circle of radius R, with R < L. In a cartesian coordinate 
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Figure 1: A crossflow fan schema. 



system x,y, the left bottom vertex of the housing is the point {—L/2, — i/2). 
The vortex waU is ideahy represented by a point, vertex of a parabola, in the 
right edge of the square, at vertical position s, with —L/2 < s < L/2 (Figl2])- 
The segment [s, L/2] is the inlet, the segment [—L/2, s] is the outlet, so that 
the fluid particles have a counter-clokwise movement into the fan. 




Figure 2: The geometry of fan. 



3 Physics and Mathematics of fan 

Now we describe the mathematical model used for simulating the air-flow into 
the fan. We consider a zero-divergence velocity field {u, v) in the 2D geometry, 
so that a scalar stream function ^ exists such that (see [6]) 



u = dy-^, V = -dx"^ 



(1) 



The expression dxV — dyU — — A'l' is called vorticity ([3])- If the case of 
a simple counter-clokwise rotating field centered at the origin, we have u = 
—•wy,v = wx (tIL), so that — A^ = 2w, where w is the rotation angular speed. 
In the case of our crossflow fan, we can also considered for fluid particles a 
negative horizontal velocity component in the upper half (0 < j/ < L/2) of the 
square and a positive one in the bottom half {—L/2 < y < 0). In the case s = 0, 
a simple way for describe this field (u/, Vf) is to consider the formulas 

Uf = —Asin{'Ky/L), Vf = (2) 

where ^ is a constant (see Figl3]). 




Figure 3: The Uf profile for L = 4.4 and A = 1. 

If s ^ 0, we use the correction Uf = —Asin[T:{s + y)/L\. Therefore, in 
the case of an estimated field {—wy + Uf,wx), the vorticity is — A^ = 2w — 
■kA/Lcos['k{s + y)/ L\. Note that, if s = and y small, using the approximation 
cosa w 1 we have 2it; — 2w cos[Tr{s + y)/ L], therefore we choose A = 1 (we 
are interested on comparison between geometries, not on absolute value of fiow 
rate) and obtain the following simple expression for the vorticity of our crossflow 
fan: 

- A* = (2w-7r/L)cos[7r(s + y)/L] (3) 

Previous formula is a differential equation of Poisson type ( j5 ) for the un- 
known function ^. On the three edges of the square that are walls for the fan 
we can impose the usual boundary conditions 

*(-L/2, y) - *(x-, -L/2) = *(x, L/2) = (4) 

because the normal component of flow along these edges is null (see [6]). 



4 Numerical resolution of Poisson equation 

If s = 0, it is possible to try a symbolic resolution of previous differential 
problem, being the equation linear and nothing that a particular solution has 
the form ^o ~ Bcos['k{s + y)/ L], with B constant. For our purposes, we apply 
a version of the numerical method of Galerkin ([8_), in the following fashion. 
We seek for a solution in the form 

N 

*(a;,2/) = X! c.n(t>ij{x,y) (5) 

where c^ are constants to be computed and (pij are the linearly independent 
functions 

^,, ={x + L/2){y + L/2)\y - L/2y (6) 

which satisfy the boundary conditions (|4]). The choice of these particular test- 
functions is justified by the fact that the flow is mainly determined by the 
position of the upper and bottom wall. Let /(y) = 2w — 'itA/Lcos['k{s + y)/ L\. 
Then the Poisson equation can be approximated as 

N 

Multiply the two members by 4>nm and integrate over the square D: 



N 



(f>n,n^4'ij) = / {f(t>nrn) (8) 
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This is a system of N'^ linear algebraic equations in the unknowns Cij which can 
be solved by usual numerical methods. 

We have computed the coefficients c^ in the case M ~ A using the methods 
built-in into the software Mathematica, ver. 7.0, for numerical integration and 
numerical resolution of linear systems. The geometrical and physical values 
which have been used for computations are impeller radius R — 2.0 cm and 
rotation rpm = 2000 {w = 27rrpm/60), L = 2{R + 0.2) cm, position of vortex 
wall pvw such that —1.5 < pvw < +1.5 cm. A graphic result of the velocity 
field in the case pvw — +1.5 is shown in Fig|4l 

5 Ecodesign optimization 

In the last years European Union has emanated some Directives about energy 
related machines, with the purpose of regulating their design and utilization 
for ecological optimizations ([2 ). In particular, some studies has revealed the 
importance of optimize the design of ventilation fans for having high energetic 
efficiency ([7]). 




Figure 4: The velocity field for pv'w=1.5 cm. 



Suppose that the rotation speed of the impeller is constant for all the possible 
values of flow rate. This is only a first quite good approximation. Then, the 
efficiency is determined by the product of the flow rate Q and the pressure head 
Ap of the fan. The pressure head is usually considered as the difference of the 
flow pressure at outlet and the flow pressure at inlet. The flow rate, defined as 
the dot product of the velocity field and the normal to the outlet surface, can 
be computed in our case with the formula 



Q- 



pvw 



L/2 



u{L/2,y)dy 



(9) 



Pressure p and velocity field v are related by Navier-Stokes equations p(Vv)v = 
— Vp, where p = 0.001 g/cm'^ is the value of air density, which we suppose to 
be constant along the flow. Once we have numerically computed v, we can 
estimate numerically the pressure along the right vertical side of the fan casing 
by a path line integral J{Vp) ■ dl (see [10]), so that if a is the ordinate of a point 
e.g. in the outlet segment [—L/2, pvw], we have 



p{a) 



dyP{L/2, y)dy 



(10) 



For our purpose, we define the pressure head of the crossflow fan has the differ- 
ence of the pressure at the middle point at outlet and the pressure at the middle 
point at inlet: 



Ap = p[l/2{pvw - L/2)] - p[l/2{pvw + 1/2)] 



(11) 



We have computed the efRciency product Ql^p for values of pvw in the interval 
[—1.5, 1.5]. Using a step of 0.1, the profile of the efficiency is shown in FiglSJ 





0.5 



1.5 



pvw 



Figure 5: Graph of the product QAp. 



The maximum is obtained for pvw = 0.6 cm, an information which were not 
easily guessed whitout a numerical analysis. From Fig|6l we can note that, with 
this geometrical configuration, the speed at outlet is almost the double of the 
speed at inlet. 




Figure 6: Profiles of u (black) and v (red) along the right side of the fan, in the 
case pvw = 0.6 cm . 
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